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Abstract 

We study the volume dependence of electric flux energies for SU(2) 
gauge theory with twisted boundary conditions. The curves interpo¬ 
late smoothly between the perturbative semiclassical results and the 
Conhnement regime. On the basis of our results, we propose that 
the Conhnement property might be caused by a class of non-dilute 
multi-instanton conhgurations. 
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1 Introduction 

Confinement has been the subject of active study since the early days of QCD. 
Wilson [1] clarihed its meaning as a possible phase of a system of Yang-Mills 
holds, and put forward the relevant order parameter, ’t Hooft added the 
characterization of other phases [2], It was clear from the start that the 
occurence or not of Conhnement was a complicated non-perturbative issue. 
Wilson’s lattice formulation [1] casted the problem into that of a peculiar type 
of statistical mechanical system, opening the way to the powerful techniques 
developed in that held of research. 

During the hrst years, most work on the subject was concentrated on 
either proving exactly the Conhnement property, or in computing by numer¬ 
ical methods the value of the string tension a ( the constant whose non-zero 
value signals Conhnement). The hrst goal has not been accomplished to date 
except in low dimensions, abelian theories [3] or in the presence of supersym¬ 
metry [4, 5]. On the numerical side, however, there is a clear signal of a 
string tension value close to the phenomenological-experimental one. 

In the last years, nonetheless, more emphasis has been set on understand¬ 
ing how the Conhnement regime arises. This amounts to having a description 
of the vacuum state of the theory and an identihcation of the relevant degrees 
of freedom that produce Conhnement. The literature is full of proposals [6], 
mostly incompatible with each other. One of the most acclaimed ideas is the 
Dual Superconductor mechanism of’t Hooft and Mandelstam [7, 8]. Lately 
this idea has been translated and studied on the lattice [9, 10]. 

Our line of attack follows the work of Luscher [11] and collaborators who 
used the volume of space as a probe for the dynamics of Yang-Mills helds. As 
one moves from small to large volumes, one is moving, without a phase tran¬ 
sition, from the region where Perturbation theory is a good approximation, 
to the one where neither the hniteness of space nor the boundary conditions 
play any role in local quantities. By observing how this process takes place 
we aim at a better understanding of the most relevant degrees of freedom in 
explaining the properties of Yang-Mills theory, most notably Conhnement. 

Our starting point differs from that of Luscher and Koller and van Baal 
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[12] in the use of twisted boundary conditions (t.b.c.). Although irrelevant 
for large volumes, the b.c.’s are crutial at small ones. Our study can be 
considered, hence, complementary to the one of the previous authors. We 
have, furthermore, reasons to believe that twisted b.c.’s are perhaps more 
efficient in showing up the most relevant conhgurations. Just on very general 
grounds, the twisted b.c.’s point towards disorder, which for a gauge system is 
Conhnement. Actually, in the large N limit t.b.c.’s seem crutial in eliminating 
the volume dependence of the string tension [13]. To exemplify what we do 
have in mind, consider the rather trivial case of the one dimensional Ising 
model. We know that this model has no ordered phase. One can explain 
this fact in terms of the kink conhgurations ( a conhguration with a single 
negative bond). The non-existence of a phase transition can be traced back 
to the non-zero probability to create a kink. Periodic b.c.’s imply an even 
number of kinks, while t.b.c.’s (antiperiodic ones) imply an odd one. The low 
temperature expansion in the hrst case starts with an ordered conhguration. 
T.b.c.’s show nonetheless the right situation from the start and the leading 
order gives us a single kink state. This pattern is followed in other less trivial 
examples. 

In this paper we are considering the case of SU(2) Yang-Mills theory 
with a spatial twist given by a magnetic hux vector m = (1,1,1) ( modulo 
2 ). An additional bonus of the twisted situation is the discreteness of the 
space of classical vacua (modulo gauge transformations), which makes the 
weak coupling expansion free of the complications of the periodic case. This 
expansion, which is a good approximation for small spatial size, has been 
studied in Refs. [14] , [15], and [16]. The hrst two references deal with 
perturbation theory alone. In the third one, we studied the role played by 
conhgurations that tunnel between the classical vacua. That work was based 
on Monte Carlo simulations in a region of small sizes, where we tested that 
the semiclassical formuli were in agreement with the data, leading to the 
numerical determination of the scale factor, which is in principle, but not in 
practice, calculable analytically. Hence, at least in the quantities that were 
studied in Ref. [16] we have a fairly good analytical understanding of the 
low Is region. 
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The present paper is to be considered a continuation of Ref. [f 6] in which 
we explore the larger values and bridge the gap towards the Conhnement 
regime. Some preliminary results were presented at the 1993 Lattice confer¬ 
ence [17]. In the next Section we give the details of the numerical method 
employed to obtain our data. In the third one we discuss the extraction of the 
physical results from the data. In the fourth section we explore some ideas 
suggested by our data, and end up presenting an scenario for understanding 
the onset of the Conhnement mechanism. In Section 5 a short summary of 
our results is presented. 

2 The Data 

In this section we present the technical details about the simulation. The 
data has been obtained on the transputer based machine RTN, located at 
the University of Zaragoza. Some results (those for smaller lattices) have 
been obtained from an 8-transputer board, similar to that contained in RTN, 
located at Universidad Autonoma. Details about RTN and the board can be 
found elsewhere [18]. The code, written in 0CCAM2 language, consists on 
a heat bath Monte Carlo simulation for SU{2) gauge helds on asymmetric 
lattices and twisted boundary conditions. Our lattices are symmetric in 
space, and the total volume is N^Nf. We have studied = 4, 6 and 8 with 
Nf = 128,128 and 64 respectively. We have been using Wilson action with 
a range of /3 values going from 2.25 up to 2.6. The specihc choices has been 
dictated by the need to cover a range of physical spatial lengths C = in 
a region where scaling is seen to apply. We have tried to obtain equal values 
of Is for different values of in order to test scaling with our data. The 
twisted boundary conditions has been introduced in our data by flipping the 
sign of the contribution of all spatial plaquettes sitting at the corner of each 
plane Vij{ni = = Ng). 

All our simulations started with several tens of thousands of sweeps in¬ 
tended to bring our statistical sample to equilibrium. The starting conhgu- 
ration has varied between the cold one, a hot (random) one and a wrongly 
cold one corresponding to no twist. This variation enables us to test thermal- 
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ization by comparing results with different initial conditions. In principle, 
the problem is more important for larger fi values, which were studied in our 
previous publication [16]. Summarising our tests, we can say that we have 
not found any signal of dependence on initial conditions in our results. 

Measurements are taken every N^kip M.C. sweeps. Typically N^kip is cho¬ 
sen on the basis of the observed autocorrelation of different quantities. Local 
quantities are quickly uncorrelated, but some of our extended Polyakov loops 
are more strongly correlated. In any case the value of N^kip which is hnally 
selected is always chosen to be larger or equal to the observed exponential 
autocorrelation times for our observables. Tests have been made that the re¬ 
sults do not change signihcantly by taking N^kip larger than the one actually 
selected to present our results. Again in this case, the stronger autocorrela¬ 
tions occur for larger values of Ng and jS. The discussion presented in Ref [16] 
is certainly sufficient to sustain our conhdence in our results. To conclude, 
we mention that the total number of M.C. iterations and measurements per¬ 
formed by our group is collected in Table 1, together with the remaining 
parameters employed in our simulation. Now we turn our attention to the 
actual measurements and observables. 

We measure ground state energies in each electric flux sector by study¬ 
ing the time dependence of correlators of the appropiate operators, namely 
Polyakov loops P. These quantities act non-trivially with respect to the 
group of singular gauge transformations in space, or equivalently carry non¬ 
zero electric flux: 

e = 03 ( 7 ) [modulo 2). ( 1 ) 

The electric flux vector e [19, 20] is a 3-vector of integers modulo 2, while 
u = (cui, a; 2 , CU 3 ) is the vector of winding numbers of the loop 7 around each of 
the 3 directions of space. In principle, one can choose any operator with the 
right quantum numbers to study these correlators. In practice, a good choice 
is required to get a reasonably good signal to background ratio. Fuzzying 
[21, 22] and Smearing [23, 24] are two techniques used in the literature to 
attain this goal. In our work we have used both techniques. The common de¬ 
nominator of both techniques is that, instead of employing straight Polyakov 
loops, one uses operators involving sums of many lattice paths with the right 
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winding. This serves two purposes: it disminishes the fluctuations, due to 
uncorrelated link fluctuations, like l/^Ah, [N^ is the number of paths in¬ 
volved), and, being a more extended operator, it is intended to increase the 
overlap with the ground state. Actually, use of these techniques is crucial to 
get reasonably small errors. 

In our version of fuzzying, we perform a renormalization group transfor¬ 
mation to a lattice of half the spatial extension. The new links being the 
normalised sum of all spatial length 2 and 4 paths between the same initial 
and final points. Straight Polyakov lines in the blocked lattice are what we 
call ’’block 1” operators. If we repeat the operation we get a lattice of one 
fourth the spatial extension and ’’block 2” operators. This can be done for 
Ns = 4: and 8. For = 6, the second blocking is built of length 3 and 5 
lines. Twisted Boundary conditions are taken into account by adding the 
different paths with appropiate signs. The sign is chosen so that all paths 
contribute equally in the large /3 limit. All of the results that we will present 
involve block-2 Polyakov loops 

Our version of the other procedure, called Smearing, is simply a local 
cooling [25, 26] step, but for each time slice independently. This procedure 
has the advantage over the previous one that it can be iterated indefinitely. 
At each ’’smearing step”, every link is replaced by the normalised sum of all 
spatial staples with the same initial and final points. Not all our simulations 
included smearing steps. Those that did, together with the number of these 
steps are shown in Table 1. 

As an alternative to smearing, some of our simulations included several 
ordinary four-dimensional cooling steps. Computing straight line Polyakov 
loops on cooled configurations can be considered as simply using some new 
cooled operators on the original configuration. Nonetheless, the new oper¬ 
ators are extended in time and the correlators of them might significantly 
change if the extension of these operators in time is comparable with their 
separation. One cooling step is, nevertheless, a local operation and hence, as 
a rule of thumb, we have limited our analysis to separations that are always 
larger than the number of cooling steps c. 

In conclusion, every Ngkip Monte Carlo iterations we measured the fol- 
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lowing quantities: , -Pi 23 ^(^)j where 6, s and c are the 

number of applied blocking, smearing and cooling steps, n is a lattice point 
and i runs over spatial directions. More precisely, we shall take 6 = 2 always, 
and cooling and smearing where never used at the same time due to memory 
limitations (See Table 1 for the characteristics of each simulation). 

We then averaged every operator of the list over spatial points. Due to 
the boundary conditions, for the |e| = 1 and 2 case, this is actually not 
projecting over an eigenstate of momentum. The Fourier descomposition of 
an operator Q{n) is 

= ( 2 ) 

S q 

where s = (e X fh) (e and m = (1,1,1) are the electric and magnetic flux 
vectors), (f is a 3-vector of integers ( modulo N^) and p = + |) is the 

momentum value. Hence, if we average over space, we get 

(3) 

S ^ s ^ i 


where h(0,(j') = and 


6 ( 1 ,^) 


i e 


■#7(3+2) 


N, sm{f-^{qP\)y 


( 4 ) 


For s = 1 all momenta appear in the sum, the smaller ones with higher 
coefficients. The only problem of using these operators instead of the pure 
eigenstates of momentum, is that the correlators may receive contributions 
from higher momentum states together with higher mass excitations. To en¬ 
sure that this does not affect our results, we have also measured the operators 
corresponding to lowest momentum, at some values of and jS. The results 
are consistent within errors, with no systematic deviation observed between 
using one type of operators or the other. 
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3 Analysis of the data 


We have measured correlators of different Polyakov loops, leading to esti¬ 
mates of the minimum energies in all electric flux sectors. In the Conhne- 
ment regime, one expects that the presence of a non-zero ’t Hooft type of 
topological electric flux, leads to the formation of a string carrying this flux. 
This physical string must have the same winding number as the Polyakov 
line which creates the electric flux. The energy of such physical string is 
given by the string tension a times the length of the string. Actually, the 
minimal length turns out to be |e| = where e^- G {0,1} is 

the component of the electric flux vector and P is the length of the torus. 
In conclusion, the minimum energies (relative to the e = 0 vacuum) in each 
sector Eg should behave as 


Eg = \e\ls(T 


(5) 


What we have actually done is to plot things in terms of the quantities 


S|e-12 



( 6 ) 


The Conhnement prediction is that, for large P, S |^2 should tend to a con¬ 
stant, equal to the string tension cr, and independent on the value of e. 
Actually, the last property (equallity of all can be interpreted as a 

consequence of the recovery of rotational invariance. 

In order to determine the energies, we measured correlations of operators 
with the appropriate quantum numbers 0{t) (the corresponding Polyakov 
loops) 

1 Nt 

C(t) = ^Y.{0{n)0{n + t)) ( 7 ) 

«=i 

from which we extracted lattice energy estimates: 


M(f+ 1) 


In 


C{t + 1) 


( 8 ) 


The smaller f, the larger the systematic errors due to contamination with 
excited states, but usually the smaller the statistical errors. For self-adjoint 
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operators one has that Mit) decreases with t. A customary compromise is to 
give as estimate of the energies the value of M[t) obtained for the smallest t 
which is compatible with larger As within errors. 

In our case, we have chosen Af(3) as an estimate for the energy. Ac¬ 
tually, even M{2) satishes the afore-mentioned criterion in most cases; and 
even when it does not, there is no clear systematics, and deviations seem to 
be due to statistical fluctuations. Nonetheless, the errors for M{2) are much 
smaller than those of Af(3), which makes the second value a safer choice. 
Additional conhdence follows from the fact that we have used different op¬ 
erators with the same quantum numbers to measure the same mass, always 
with consistent results. The different operators are the ones associated with 
different smearing steps and both the momentum eigenstates, as the space 
averages, when available. When cooling was used, according to the reasoning 
of the last section, we have actually taken Af(3 + c) as our estimate. Again 
compatibility within errors was attained. Our conclusion is that, either there 
is little contamination of higher states, or else the latter have energies that 
differ from the ground state ones by less than the errors. 

In order to give the reader an idea of the quality of our data, we have 
collected in Table 2 a set of determinations for several operators and values 
of Ws, fi and time separations t. Errors have been estimated by partitioning. 
We show values of Sj- in physical units ( see below) for small ( 0.671 fm.), 
intermediate ( 0.910 fm.) and large volumes (1.233 fm.). The data at /3 = 2.3 
and Ns = 4: represents two independent runs with up to 3 cooling steps, but 
differing by the use of momentum eigenstate operators ( MEO) or simple 
averages. The intermediate volume result given in Table 2 employed up 
to 9 smearing steps. Einally, the large volume one again consisted on two 
independent runs, one with cooling and the other with smearing. The c = 
s = 0 Mixed result combines the statistics of both runs. All our results are 
represented by converting energies into physical units by use of the empirical 
formula 

a{p) = 400 f^- (9) 

which applies in our region of /3 values according to the results of Ref. [27]. 
The absolute normalization is hxed by convention to cr = 5/m“^. 
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In Fig. 1, E |^|2 is converted into physical units and plotted as a 

function of Is- We have added the results of Ref. [16] , which affect only 
e = (1,1,1) and < 0.7 fm. The plotted values follow from the Af(3 + c) 
estimates of the energies. The best choice for the number of coolings c or 
the number of smearing steps s, is selected. Notice that in some cases data 
from different lattice sizes give very similar values of Is- Our data shows that 
scaling is satished to within a few percent. 

At low values of the predictions of the semiclassical approximation are 
also shown. Both Si and S 2 receive a contribution from perturbation theory, 
going like The coefficients were computed in Ref. [15]. S 3 is zero to all 
orders in perturbation theory, and becomes non-zero as a result of tunneling 
mediated by the Q = ^ instanton [28]. In Ref. [16] we showed that indeed 
the data follow the predictions of the semiclassical dilute gas approximation. 
This determines the dependence of S3 up to a multiplicative prefactor 
AslWs), which, in the absence of an analytical computation, has been htted 
to our data. In the case of Si and S 2 , we also allowed for a non-perturbative 
piece of the same shape but with a free prefactor. 

Hence, the curves shown in Fig. 1 correspond to the functions: 

W(iV„ /3) = + B,{N,) + ^ , (10) 

where S{Ns) is the lattice action of the Q = ^ instanton, Bi{Ns) and Ci{Ns) 
(i = l,2)are given in Ref. [15], and is htted to our data. Actually, 

for |ep ^ 3, the dominant term is the one in the middle, where Bi = B 2 
is equal to 4.14, 4.30 and 4.36 for = 4, 6 and 8 respectively. The next 
perturbative piece has the right tendency for the splitting between Si and S 2 
shown by the data, but its magnitude is very small ( 1 % ). Finally, the non- 
perturbative piece gives contributions ranging from 3 to 10 % as one moves 
from Is = 0.4 to 0.7fm. Due to the size of the errors, setting A 2 = A 3 = 0 
still gives rise to a reasonably good ht, and the curves sit in between the Si 
and S 2 data. The many curves seen in Fig 1. come from the different values 
of Ng. Their closeness is again a manifestation of scaling for our data. 

In summary, we can say that the semiclassical predictions describe our 
data fairly well, conhrming that we have a reasonably good understanding 
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of the dynamics occuring for volumes smaller than 0.7fm, which is the bor¬ 
der of the dilute gas approximation region. On the opposite side, for large 
/s, Conhnement tells us that the different E’s should merge into a unique 
horizontal line of height 5/m“^. Indeed our data show this tendency, and 
furthermore the interpolation between both regions is fairly smooth. This 
behaviour is the main result of our data, and in the next section we will try 
to understand how this could take place. 

4 An scenario for Confinement 

The purpose of this section is to analyse the implications of our Monte Carlo 
results on the nature of the SU{2) gauge theory vacuum and Conhnement. 
We expect that these hndings could be translated to SU{3) theory and QCD. 
The most important conclusions from our data are two. First, it determines 
the scale at which the hnite volume theory approaches the inhnite volume 
one to within a 20% accuracy. Our quantities S |^|2 actually point towards 
1.0 — 1.2 fm as the appropiate scale. Furthermore, precisely at this scale 
rotational invariance is recovered with the same level of accuracy. Second, 
the approach to the inhnite volume value is smooth. There seems to be no 
scale where a qualitative change of behaviour is observed. This is reinforced 
if we take into account that also the glueball spectrum is qualitatively alike 
for large and small volumes [29]. This is in contrast with the situation for 
periodic boundary conditions [30]. 

Hence, it is natural to expect as well a smooth transition between the 
conhgurations which dominate the path integral for different values of Ig. Of 
course, it is logically possible that a qualitative change in the structure of 
vacuum does not lead to drastic changes in the values of S. We nevertheless 
would like to explore the implications of our more natural hypothesis. In¬ 
deed, we have a very good understanding of the structure of the vacuum for 
Is < 0.7 fm. It consist on a gas of dilute twisted periodic instantons [28] 
modulated by the gaussian (and higher) perturbative corrections. This is 
seen with textbook clarity for very small and describes perfectly well our 
data for Ig < 0.7 fm as seen in Fig. 1 and Ref [16]. Hence, we are lead to 
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conjecture that the descendants of the twisted periodic instantons do popu¬ 
late the Yang-Mills vacuum. This is very shocking at hrst sight, and it is the 
purpose of this section, to clarify how this can happen and to extract some 
physical consequences from this fact. We will end up proposing an scenario 
which has several appealing features. One of these, is the possibility of ex¬ 
plaining both the value of the string tension and the topological susceptibility 
in terms of the same classical conhgurations. Even more, it is possible that 
chiral symmetry breaking also originates from them [31]. 

Before proceeding to describe what can happen for ~ 1.2 /m, let us 
consider again the small region. The number of Q = | instantons per unit 
time n/, keeps increasing with Is- Indeed, the time extent of each instanton 
is /s, so that the dilute gas approximation breaks down when 

1 / ^ 

nj ~ . (11) 


Given the relationship between S 3 and n/, we have that this occurs at 
such that 

Mh ) ^ -^=2 , ( 12 ) 

and this can be seen to hold for ~ 0.7 fm. Hence, beyond this point 
the instantons are very close to each other and one has to take into account 
their interaction. Within the dilute gas approximation there is a relationship 
between the topological susceptibility and S 3 which is as follows 


= ^S 

^ 8/2 ■ 


(13) 


Thus, at the border of the dilute gas approximation, one has y (fm) ^ 
which is in agreement with what other groups have measured in their Monte 
Carlo simulations [32]. 

What happens beyond 0.7/m ? There are several possibilities. Q = ^ 
instantons might fuse and make ordinary Q = 1 instantons. However, the 
latter are ineffective in producing Conhnement and hence, if this was the 
case we should observe a decrease in the string tension precursor S3. It is 
very unnatural to think that some other mechanism could come in precisely 
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at that point and keep the things slowly varying. There is indeed something 
else which we know must happen: some other conhgurations can start to play 
a role in the path integral. In particular one can consider the descendants 
of our twisted periodic instanton. These conhgurations are periodic in space 
but with periods which are fractions of Is- Let us note this period by d. We 
then have that l^/d is an integer . Notice that these conhgurations are special 
types of multi-instanton conhgurations on the torus ( or in oo space). They 
exist both for twisted or periodic boundary conditions, with only an overall 
constraint on their number distinguising different b.c.’s. 

How special are these conhgurations? If we have one such periodic con- 
hguration with N lumps, it has a total topological charge of Q = N/2. 
According to the index theorem [33] there are 8 Q = AN linearly indepen¬ 
dent deformations of these solutions into other multi-instanton solutions. 
Only 4 of the total AN maintain the periodic arrangement. Notice that the 
number of degrees of freedom of the multi-instanton space equals that of a 
4-dimensional gas of particles. We propose to view the multi-instanton space 
precisely as a gas of Q = ^ lumps, where one needs only to give the centers 
of these lumps to specify the conhguration uniquely. These lumps behave 
as elastic balls or balloons, whose size and shape adjusts appropiately given 
their location. The lumps are not isolated or dilute ( in all directions); they 
cannot be, since their size grows with separation. It is rather a liquid than 
a gas. This description is obviously possible and valid in the neighbourhood 
of the periodic conhgurations. For large deformations it must fail. When 
two lumps are pushed into each other, they fuse into a single Q = f lump, 
and their 8 parameter space now gives the position, size and orientation of 
the lump (as for ordinary instantons). By generating conhgurations on the 
lattice, we have actually observed how there is a continuous path from the 
two Q = ^ lump structures to the Q = f ones [34]. 

Summarising, we are proposing to consider a hnite measure part of the 
multi-instanton space, which can be described as a liquid of Q = ^ lumps. 
We will label such space by the mean-separation d among lumps, which is 
limited by the hnite spatial size and the boundary conditions. We will de¬ 
duce the properties of these conhgurations from those of the periodic ones. 
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of which they are deformations. An analytic and numerical understanding of 
this parametrization of multi-instanton space and its properties is a purely 
classical problem. By now only numerical evidence in favour has accumu¬ 
lated, but the problem is currently under study. 

How do these conhgurations become relevant at larger volumes? To study 
this point, we saturate the path integral with this kind of conhgurations. The 
structures are naturally 3 dimensional for small volumes, and so we will make 
the counting in a way which privileges time. We have 


n{d) 


^-n(d)F(d,ls) 


( 14 ) 


where n{d) is the number of space-periodic conhgurations with period d 
(a fraction of Ig). The relative weight of one such conhguration is 
If we make the simplifying hypothesis (which we will comment upon later) 
that 

^-F{dM) ^ ^^-F{d,d)yif r ^ 

this relative weight is a function of exp{—F[d^ d)} which can be extracted 
from our data in the region where the dilute gas approximation applies: 





( 16 ) 


with a ~ 2 and Ig — 0.6 — 0.7 fm. This formula can be assumed to be valid, 
by extrapolation, for Ig slightly larger than Ig as well. In the thermodynamic 
limit [Ig —)• oo with d hxed), the relative weight peaks at d = d given by 
the maximum of the function 


G(d) = 


( 17 ) 


Indeed, the previous function has a maximum which is independent of a, 
given hy d = IgC^d > Although the precise relation between d and Ig 
depends on the details of the computation, the existence of a maximum and 
the fact that d > Ig, can be easily understood. For d < Ig it is indeed 
more probable to have the perturbative vacuum than to have an instanton 
of size d in one cell of size d'^. Actually, since these instantons are periodic 
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in space, the larger the volume the more suppressed such a conhguration 
is. The tendency changes when d > Is- Then the more instantons the more 
probable. However, the probability of one instanton increases with d, but 
the number of instantons that ht into a hxed spatial volume decreases with 
d. Hence there is a maximum within this region. 

Although our approximations are presumably fairly crude, the conclusion 

is quite robust to modihcations: In the the thermodynamic limit (/^ -)• oo), 

only one of the family of periodic conhgurations with period d dominates the 
path integral, namely that with d = d the maximum of G{d). This leads 
to the following description of what could happen as the volume of space 
increases from 0 to oo. For /^ < d ~ the dominating conhguration of in¬ 
stantons has d = Is and we have a dilute gas of these conhgurations. Beyond 
Is ^ d we start having conhgurations with periods d = Is/n (n integer). The 
most probable one actually corresponds to d = d independently of Is- This 
explains how one is hnally going to obtain results which are volume indepen¬ 
dent for large volumes. The results will also be independent on the boundary 
conditions imposed in space: periodic or twisted. Actually, boundary condi¬ 
tions only hx the evenness or oddness of the number of periodic instantons 
contained in space, which is irrelevant for large numbers. 

Now we have to face the meaning and justihcation of formula (15). Ac¬ 
tually, T(d, Is) gives the free energy relative to the vacuum of the periodic 
instanton conhguration. In the semiclassical approximation one has: 

( 18 ) 

where the exponential term is exp{ — ^} and Sd is the classical action of the 
conhguration. In our case this action is dvr^ for each cell of size d^, in agree¬ 
ment with (15). Hence, we simply need to prove and analogous property 
for the prefactor A{d^ls). This quantity incorporates the integration over 
gaussian huctuations around the classical conhguration and the integration 
over the collective coordinates of the zero-modes. For gaussian huctuations 
of high frequency (;> 1/d) it is reasonable to expect that the huctuations 
around each conhguration are insensitive to the environment in which they 
are placed. This translates again into a dependence of type (15). Deviations 
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from this formula are to expected for frequencies G and our hypothesis 

amounts to assuming that their effect does not disturb expression (15) too 
much. Concerning zero-modes, notice that to compute A{d^ d) one has to 
consider the only four zero-modes associated to space-time translation of the 
conhguration. To justify (15) we should have 4 zero-modes for each cell of 
size which is precisely what we expect from the index theorem. Notice, 
that irrespective of C, one expects that each collective coordinate associated 
to a zero-mode is of range d. We stress that, since we are integrating over 
the zero-mode space, we are actually considering not-only spatially periodic 
conhgurations, but also the other non-periodic conhgurations: The whole set 
of positions of the gas of lumps in our scenario. Our main assumption in this 
case is that gaussian fluctuations around these non-periodic conhgurations 
are similar to those for periodic ones. In conclusion, expression (15) is the 
most reasonable approximation to be made in the absence of an exact ex¬ 
pression. We expect that this assumption gives qualitatively correct results 
although quantitative corrections are likely. 

The previous considerations imply that the main free parameter in our 
instanton liquid model of the vacuum, the mean separation d, is hxed dy¬ 
namically in the average to a value close to 0.7/m. We will see that such 
a class of conhgurations gives rise to properties of the vacuum which are 
in reasonable qualitative agreement with the known properties of the SU(2) 
gauge vacuum. First we will consider the topological susceptibility y and 
next the string tension. 

Assuming as usual that there are equal probabilities for Q = ^ instantons 
and anti-instantons, the topological susceptibility turns out to be equal to 

X = ^ - l/m“^ . (19) 

This value is in agreement with the results of Monte Carlo measurements 
of other authors [32]. Actually, since our instantons are non-dilute, there is 
no reason why the probability of instantons and anti-instantons should be 
independent. There seems to be some evidence in the direction of a higher 
probability of alike charges to cluster [35], but this might be also due to 
the distortion caused by cooling. In any case, at our qualitative level, a 
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reasonable amount of clustering can be compensated with a small change in 

d. 

Finally, we should touch upon a crucial point. Does the proposed liquid 
give rise to a non-vanishing string tension? Although it might be possible to 
show such a behaviour by numerical methods, we do not yet have results in 
this direction. Here, we will try to make a crude analytical approximation 
which is expected to give the right order of magnitude. Again our consider¬ 
ations will make use of our knowledge of the periodic Q = ^ conhguration in 
a 4 dimensional twisted box of size d: our prototype of an individual lump 
[36]. If we measure the square Wilson loop around the boundaries of the box 
centered at the lump, we hud that if the position of the loop in the trans¬ 
verse directions coincides with the center of the loop, then the value of the 
Wilson loop is —I. On the other hand if the Wilson loop lies in between two 
lumps, the result is /. At intermediate positions the loop gives an element of 
the group interpolating between both, but not situated along a single U{1) 
subgroup of SU(2). Since computing Wilson loops in classical conhgurations 
is notably difficult, we will make use of the following customary abelianizing 
approximation. We will assume that the loop is only either —I or /. The hrst 
case occuring when the loop passes within a distance xd from the center of 
the lump. Otherwise the result is /. Typically this approximation amounts 
to replacing a continuous function going from -1 to 1, by a step function. 
The parameter x should lie between 0 and |. Now, if we have a large Wilson 
loop, its value in a given conhguration is equal to ( — 1)^, where M is the 
number of lumps which are within a distance xd from the minimal surface of 
the loop. In our gas picture, M has a Poissonian-binomial distribution. The 
value one gets for the string tension is: 

2 

^ ( 20 ) 

where g{x) depends on the details of the distribution: for Poissonian, 

— ln[l — 8a:^)/2 for a binomial one. Notice that for g(x) of the order of 1, one 
gets a value of the string tension which is close to the one measured for large 
lattices. Since the string tension depends on d , a relatively small change 
in d can adjust for corrections to the approximation. 
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5 Conclusions 

In this paper we have presented a follow up of the evolution of the value of 
electrix flux ground state energies as a function of the spatial volume. Our 
data nicely show how and where the predictions of Conhnement are obtained 
for large lattices. No intermediate scales are appreciated in our data and at 
sizes of /s ~ 1 — 1.5/m both rotational invariance and Conhnement pre¬ 
dictions are well reproduced. Since an analytic understanding of the low 
region is known [16], we are led to speculate on the way that the correspond¬ 
ing structure of the vacuum smoothly evolves into the conhning vacuum. 
An appealing scenario arises in which an important fraction of the space of 
multi-instanton conhgurations can be described as a liquid of (non-isolated) 
(5 = 1 instantons. From results of our previous work [16], we argue that the 
average density of lumps in the vacuum is dynamically hxed to be (0.7/m)“‘^ 
approximately. With standard assumptions, we show that this gas could lead 
to values of the topological susceptibility y, and of the string tension a which 
are close to the ones known to hold for SU(2) Yang-Mills theory. 

Our scenario has in some respects striking similarities with the one pro¬ 
posed by Callan, Dashen and Gross [37]. First of all, the main constituents 
of our vacuum have | topological charge, and can be combined in pairs and 
smoothly deformed to make ordinary instantons. These properties are shared 
by the corresponding objects in the case of Ref. [37], which are called merons 
[38]. Analogously, the mechanism by means of which merons produce Con¬ 
hnement is virtually the same that we have employed in computing the string 
tension for our liquid. Nonetheless, there are also important differences be¬ 
tween our Q = \ instantons and merons. In our case, the conhgurations 
are smooth everywhere, while merons are singular solutions. It might be 
possible nonetheless, that one could establish some connection between our 
conhgurations and merons. 

Unfortunately, unlike the meron case, there is no analytic formula for the 
Q = \ instantons, although they are quite well known numerically. This 
problem difhcults progress on the analytic side. It is nevertheless possible 
to study these problems numerically on the lattice. The Q = \ lumps are 
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seen to be present in our Monte Carlo generated configurations [39]. We are 
currently investigating the structure of the cooled SU(2) gauge vacuum in 
order to prove or disprove this or other models of Conhnement. Some of our 
results are strikingly consistent with our model [40]. 
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Table captions 

1. Table 1: 

We show the total list of simulations that we have performed, indi¬ 
cating the value of fi and the initial conhguration, the number of 
measurements, the number of Monte Carlo sweeps between two suces- 
sive measurements N^kip and the number of thermalization sweeps Nq. 
We indicate how many cooling steps and smearing steps, if any, we 
have performed and whether we measure the three fluxes or only the 
one with |ep = 3. 

2. Table 2: We give the values of S |^2 and errors for some of our simulation 
points [Ng and fi). The columns labelled indicate that the result is 
extracted from correlations at distances c + n and c + n + 1 , where c is 
the number of coolings. For each simulation, different estimates of S 
from different operators is shown (s stands for smearing step) . 


Figure captions 

f. Figure 1: 

The values of S |^2 (dehned in Fq. ( 6 )) are plotted as a function of C. 
The different symbols used for different values of are shown in the 
plot. The semiclassical curves are described in the text and from top 
to bottom correspond to Si = 6 and 4, S 2 iV^ = 6 and 4 and S 3 
Ws = 4, 6 and 8 . 
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Table 1: 


Ns 

/3 

Genesis 

Meas. 

^skip 

iVo 

Type 

Flux 

8 

2.325 

cold 

568 

500 

40000 

s=0,3,6,9 

1,2,3 

8 

2.325 

cold 

1336 

250 

324000 

c=0,l,2,3 

1,2,3 

8 

2.360 

cold 

750 

500 

40000 

s=0,3,6,9 

1,2,3 

8 

2.385 

hot 

310 

500 

40000 

s=0,3,6,9 

1,2,3 

8 

2.415 

cold 

300 

500 

40000 

s=0,3,6,9 

1,2,3 

6 

2.330 

cold 

644 

200 

40000 

s=0,3,6,9 

1,2,3 

6 

2.369 

cold 

1211 

100 

56000 

s=0,3,6,9 

1,2,3 

6 

2.395 

hot 

740 

200 

40000 

s=0,3,6,9 

1,2,3 

6 

2.425 

cold 

1127 

100 

40000 

c=0,l,2,3 

3 

6 

2.460 

cold 

620 

250 

40000 

c=0,l,2,3 

1,2,3 

6 

2.483 

hot 

335 

500 

42500 

s=0,3,6,9 

1,2,3 

4 

2.300 

hot 

818 

100 

40000 

c=0,l,2,3 

1,2,3 

4 

2.300 

hot 

541 

100 

121800 

c=0,l,2,3 

1,2,3 

4 

2.300 

cold 

3920 

50 

20000 

c=0,l,2 

3 

4 

2.340 

hot 


100 

70000 

s=0,3,6,9 

1,2,3 

4 

2.340 

cold 

3920 

50 

20000 

c=0,l,2 

3 

4 

2.360 

cold 

686 

100 

120000 

s=0,3,6,9 

1,2,3 

4 

2.380 

hot 

1692 

50 

300000 

c=0,l,2,3 

1,2,3 

4 

2.380 

hot 

642 

50 

384600 

c=0,l,2,3 

1,2,3 

4 

2.400 

hot 

1692 

100 

40000 

c=0,l,2,3 

1,2,3 

4 

2.425 

hot 

484 

100 

40000 

c=0,l,2,3 

1,2,3 

4 

2.440 

hot 

319 

250 

20000 

c=0,l,2,3 

1,2,3 

4 

2.475 

hot 

219 

500 

40000 

s=0,3,6,9 

1,2,3 
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Table 2: 


Ns 

/3 

Type 

|e|2 = 1 

|e|2 = 1 

|e|2 = 2 

|e|2 = 2 

|eT = 3 

|eT = 3 




h 

^3 

h 

^3 

h 

^3 

8 

2.325 

s=0 

5.20(13) 

4.80(31) 

5.84(29) 

4.9(1.5) 

4.77(61) 

- 

8 

(1.233) 

s=3 

5.27(16) 

4.83(43) 

5.18(28) 

5.1(1.3) 

5.10(61) 

13 ( 11 ) 

8 


s=6 

5.42(18) 

4.59(42) 

5.20(32) 

5.2(1.9) 

5.04(67) 

- 

8 


s=9 

5.47(21) 

4.66(50) 

5.36(37) 

5.2(1.8) 

5.10(70) 

- 

8 

2.325 

c=0 

5.029(73) 

4.62(22) 

5.30(19) 

4.28(55) 

4.43(33) 

3.9(2.0) 

8 

(1.233) 

C=1 

4.83(12) 

4.75(39) 

4.59(27) 

3.24(96) 

3.61(46) 

3.3(1.8) 

8 


c=2 

4.79(26) 

4.8(1.1) 

3.83(74) 

2 .9(1.7) 

2 .8(1.1) 

3.5(1.8) 

8 


c=3 

4.85(78) 

- 

2 .7(1.5) 

- 

3.8(2.3) 

- 

8 

MIXED 

c=0 

5.098(63) 

4.75(17) 

5.36(13) 

4.46(44) 

4.52(21) 

5.01(30) 

8 

2.415 

s=0 

6.36(15) 

6.41(28) 

6.64(20) 

6.40(53) 

4.18(20) 

3.69(48) 

8 

(0.910) 

s=3 

6.07(14) 

6.64(37) 

6.48(17) 

5.63(45) 

3.95(17) 

3.44(36) 

8 


s=6 

6.08(17) 

6.62(40) 

6.36(18) 

5.75(49) 

3.69(20) 

3.56(28) 

8 


s=9 

6.20(18) 

6.62(44) 

6.28(23) 

5.74(61) 

3.57(14) 

3.73(34) 

4 

2.300 

c=0 

9.25(18) 

9.90(58) 

8.40(25) 

7.85(62) 

2.837(51) 

2.821(77) 

4 

(0.671) 

C=1 

8.88(40) 

8.26(88) 

7.93(50) 

7.3(1.3) 

2.791(72) 

2 .88(12) 

4 


c=2 

8.35(71) 

8 .7(1.6) 

7.2(1.1) 

4.6(2.2) 

2.87(10) 

2.99(26) 

4 


c=3 

8 .6(1.3) 

8.4(3.1) 

4.4(1.9) 

- 

2.93(26) 

2.96(44) 

4 

2.300 

c=0 

9.03(27) 

8.58(62) 

7.92(31) 

7.9(1.1) 

2.847(62) 

2.667(77) 

4 

MEO 

C=1 

8.80(44) 

9.5(1.2) 

8.10(75) 

7.2(2.5) 

2.765(67) 

2.724(98) 

4 


c=2 

9.15(97) 

11.9(3.1) 

6 .9(1.6) 

4.7(5.7) 

2.780(92) 

2.82(18) 

4 


c=3 

10.7(2.0) 

13.3(8.9) 

6 .9(7.5) 

- 

2.81(18) 

2.87(36) 


25 
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